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ABSTRACT 



Aims. We investigate observational constraints on the generalized Chaplygin gas (GCG) model including the gamma-ray bursts 
(GRBs) at high redshift obtained directly from the Union2 type la supernovae (SNe la) set. 

Methods. By using the Markov Chain Monte Carlo method, we constrain the GCG model with the cosmology-independent GRBs, 
as well as the Union2 set, the cosmic microwave background (CMB) observation from the Wilkinson Microwave Anisotropy Probe 
(WMAP7) result, and the baryonic acoustic oscillation (BAO) observation from the spectroscopic Sloan Digital Sky Survey (SDSS) 
data release 7 (DR7) galaxy sample. 

Results. The best-fit values of the GCG model parameters are A s =0.74751g;°|«(la-)!°;^(2tr), or=-0.0256+°;}™(lo-)^™(2o-), and 
the effective matter density Q. m = 0.2629^^ Q]5|(lcr)^Q Q|||(2cr), which are more stringent than previous results for constraining GCG 
model parameters. 
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1. Introduction 

The original Chaplygin gas (CG, Kamenshchik et al. 2001) and 
generalized Chaplygin gas models (GCG, Bento et al. 2002) 
have been proposed as possible explanations of the acceleration 
of the current universe, with the equation of state as follows 



Pgcg - — 



Pgcc 



(1) 



where A and a are two parameters to be determined. For the 
case a = 1, it corresponds to the original Chaplygin gas 
(Kamenshchik et al. 2001); if a = 0, it acts as the cosmolog- 
ical constant (A). Considering the relativistic energy conserva- 
tion equation in the framework of Friedmann-Robertson- Walker 
(FRW) metric, we can obtain 



Pgcg - Pgcg 



+ (1 - A s )a 



-3(1 +a) 



r 



(2) 



where A s = A/p^ CQ , pgcg.o is the energy densities of the 
GCG at present, and the scale factor is related to the redshift 
by a = 1/(1 + z). From Equation (2), the striking property of the 
GCG can be found that the energy density behaves as dust like 
matter at early times; while it behaves like a cosmological con- 
stant at late times. Therefore, the GCG model can be regarded as 
a derivative of the unified dark matter/energy (UDME) scenario 
(Bento et al. 2004). Until now, the GCG model has been con- 
strained using many different types of observational data, such 
as Type la supernovae (SNe la) (Fabris et al. 2002; Makler et al. 
2003a; Colistete et al. 2003; Silva and Bertolami 2003; Cunha 
et al. 2004; Bertolami et al. 2004; Bento et al. 2006; Wu and 
Yu 2007a), cosmic microwave background (CMB) anisotropy 
(Bento et al. 2003a, 2003b; Bean and Dore 2003; Amendola et 



al. 2003), the angular size of the compact radio sources (Zhu 
2004), the X-ray gas mass fraction of clusters (Cunha et al. 
2004; Makler et al. 2003b), the Hubble parameter versus redshift 
data (Wu and Yu 2007b), large-scale structure (Bilic et al. 2002; 
Multamaki et al. 2004), gravitational lensing surveys (Dev et al. 
2003, 2004; Chen 2003a, 2003b), age measurements of high-z 
objects (Alcaniz et al. 2003) and lookback time of galaxy clus- 
ters (Li, Wu and Yu 2009); as well as various combinations of 
data (Wu and Yu 2007c; Davis et al. 2007; Li, Li and Zhang 
2010; Xu and Lu 2010). 

Gamma-ray bursts (GRBs) have been proposed as dis- 
tance indicators and regarded as a complementary cosmologi- 
cal probe of the universe at high redshift (Schaefer 2003; Dai 
et al. 2004; Ghirlanda et al. 2004; Firmani et al. 2005, 2006; 
Liang and Zhang 2005; Ghirlanda et al. 2006; Schaefer 2007; 
Wang et al. 2007; Wright 2007; Amati 2008; Basilakos and 
Perivolaropoulos 2008; Mosquera Cuesta et al. 2008a, 2008b; 
Daly et al. 2008). Owing to the lack of a low-redshift sample, 
the empirical luminosity relations of GRBs had been usually 
calibrated by assuming a certain cosmological model with par- 
ticular model parameters. Liang et al. (2008) presented a com- 
pletely cosmology-independent method to calibrate GRB lumi- 
nosity relations with the luminosity distances of GRBs at low 
redshift interpolated directly from SNe la or by other similar ap- 
proaches (Liang and Zhang 2008; Kodama et al. 2008; Cardone 
et al. 2009; Gao et al. 2010; Capozziello and Izzo 2010). 
Following the cosmology-independent calibration method, the 
derived GRB data at high redshift can be used to constrain cos- 
mological models by using the standard Hubble diagram method 
(Capozziello and Izzo 2008; Izzo et al. 2009; Wei and Zhang 
2009; Wei 2009; Qi et al. 2009; Wang et al. 2009a, 2009b; Liang, 
Wu and Zhang 2010; Wang and Liang 2010; Liang, Wu and Zhu 



2 



N. Liang, L. Xu and Z.-H. Zhu: Constraints on GCG model including GRBs 



2010; Wei 2010a, 2010b; Liang and Zhu 2010; Demianski et al. 
2010). Bertolami and Silva (2006) first studied the GCG model 
by considering the use of GRBs at 1.5 < z < 5 calibrated with 
the bursts at z < 1.5 as distance markers. The joint analysis with 
the GCG model of the cosmology-independent GRB data set ob- 
tained in Liang et al. (2008) can be found in Wang et al. (2009a) 
and Freitasa et al. (2010). 

Liang, Wu and Zhu (2010) calibrate GRBs data at high red- 
shift directly from the Union2 compilation of 557 SNe la data 
set (Amanullah et al. 2010) and constrain the Cardassian mod- 
els by combining the updated GRB data with the joint obser- 
vations. In this paper, we use the Markov Chain Monte Carlo 
(MCMC) technique to constrain the GCG model from the latest 
observational data including the updated distance moduli of the 
GRBs at high redshift obtained directly from the Union2 set. We 
combine the GRB data with both the joint observations, such 
as the Union2 set, the CMB observation from the Wilkinson 
Microwave Anisotropy Probe (WMAP7; Komatsu et al. 2010) 
result, and the baryonic acoustic oscillation (BAO) observation 
from the spectroscopic Sloan Digital Sky Survey (SDSS) data 
release 7 (DR7) galaxy sample (Percival et al. 2010). 



data as far as cosmological constraints go. The acoustic scale can 
be expressed as 



la =n 



£X: ll2 sinn[£ll> 2 £* ^]/H 



(6) 



where r s (z„) = Hq~ 1 (°° c s (z)/E(z)dz is the comoving sound 
horizon at photo-decoupling epoch. The shift parameters can be 
expressed as 

dz 



R 



n 1/2 tr 1/2 



,1/2 



f 

Jo 



(7) 



recombination can be given by (Hu and 



The redshift of 
Sugiyama 1996) 

z, = 1048[1 +0.00124(Q/,/z 2 r a738 (l 

2\-0.238/ 



where gl = 0.0783(Q/,/z 2 ) 

!,2s.l.8h-l 



• 8l (Cl MO h z ) g2 )], 

r,2\-0.763\-l 
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5 (1 +39.5(£Wr a/M r 1 andg 2 = 
0.560(1 + 21 .l(CVz z ) 1M r'- From the WMAP7 measurement, 



the best-fit values of the data set (/„, R, z„) are (Komatsu et al. 
2010) 



CMB 



2. Observational Data Analysis 
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The Union2 compilation consists of data for 557 SNe la 
(Amanullah et al. 2010), and we use the 69 GRBs data com- 
piled by Schaefer (2007). Following Liang, Wu and Zhu (2010), 
we use the updated distance moduli of the 42 GRBs at z > 1.4, 
which calibrated with the sample at z < 1.4 by using the lin- 
ear interpolation method from the Union2 set. For more details 
for the calculations for GRBs, we refer to Liang et al. (2008^ CMB " 1 - 
and Liang, Wu and Zhang (2010). Constraints from SNe la and 
GRB data can be obtained by fitting the distance moduli /i(z). A 
distance modulus can be calculated as 



The x 2 value of the CMB observation can be expressed as 
(Komatsu et al. 2010) 

Cr 



AP, 



CMB 



-CMB 'APcMB, 



(10) 



where AP CM b = Pcmb 
covariance matrix is 

2.305 29.698 -1.333 
29.698 6825.270 -113.180 
-1.333 -113.180 3.414 



Pcmb, and the corresponding inverse 



(11) 



fi = 5 log — — + 25 = 5 log 10 D L - ju , 
Mpc 



(3) 



For the BAO observation, we use the measurement of the 
BAO distance ratio (d z ) at z = 0.2 and z = 0.35 (Percival et al. 
2010), which can be expressed as 

r s (zd) 



d 7 = 



(12) 



where fi () = 5 log 10 [// /(100km/s/Mpc)] + 42.38, and the lumi- 
nosity distance D L can be calculated using 



D L = H d L = (1 + z)£r 1/2 sinn|Q, 



,1/2 



I 



dz' 



E(z') 



Dv(zbao) ' 

where the distance scale Dy is given by (Eisenstein et al. 2005) 

ll/3 



(4) 



(13) 



where sinn(x) is sinh for Q.^ > 0, sin for Qk < 0, and x for 
Ok = 0, and E(z) = H/Hq, which is determined by the choice of 
the specific cosmological model. The x 2 value of the observed 
distance moduli can be calculated by 



Dv(ZBAo) = io[^(I wfi 

and r s (zd) is the comoving sound horizon at the drag epoch at 
which baryons were released from photons, Zd can be given by 
(Eisenstein & Hu 1998) 



Zd = 



1291(0 M o/! 2 ) 



■[(l + ^cv* 2 / 2 )], 



[1 + 0.659(Qmo^ 2 ) 



0.828 



(14) 



2 = y [jUobsfc) -j»fa)] 2 
i=l U M,i 



(5) 



BAO 



where fi bs(zd are the observed distance modulus for the SNe 
la and/or GRBs at redshift z, with its error cr^.; yu(z,) arep 
the theoretical value of the distance modulus from cosmo- 
logical models. Following the effective approach (Nesseris 
and Perivolaropoulos 2005), we marginalize the nuisance 
parameter /io by minimizing xlt - C - B 2 /A, where 
A = 1,1 /a%, B = S^obsfaO-Slog^Di]/^ 2 , an d C = 
Z[^ob s (Zi)-5log 10 D L ] 2 /crl. 

For the CMB observation, we use the data set including the 
acoustic scale (l a ), the shift parameter (R), and the redshift ot BAO _1 = 
recombination (z„), which provide an efficient summary of CMB 



where b\ = 0.313(£Wr) y [l + 0.607(Q M0 ^) T and 
b 2 = 0.238(fWz 2 ) - 223 - From SDSS data release 7 (DR7) 
galaxy sample, the best-fit values of the data set (do.2, ^0.35) are 
(Percival et al. 2010) 

Jo 2 \ _/ 0.1905 ±0.0061 \ 
do.35 J ~ ^ 0.1097 ±0.0036 J' 

The x 2 value of the BAO observation from SDSS DR7 can be 
expressed as (Percival et al. 2010) 



(15) 



2 

^BAO 



_ APg AQ CBAO 'APbao, 



where the corresponding inverse covariance matrix is 

30124 -17227 \ 
-17227 86977 ' 



(16) 



(17) 
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3. CONSTRAINTS ON THE GCG MODEL VIA MCMC 
METHOD 

We consider a flat universe filled with the GCG component and 
the baryon matter component. From the Friedmann equation 
H 2 = (87rG/3)(p b + Pgcg), we find that 

H 2 

E 2 (z;A s ,a)= — = £2 b (l + z) 3 + (1 - Q b ) 
H o 

x[A s + (I - A s )(l + zf (1+a) ]~ , (18) 

where Q b represents the fractional contribution of baryon matter. 
The effective matter density in the GCG model can be given by 
(Bento et al. 2004; Wu and Yu 2007c) 

O m = Qb + (l-Ob)(l-A,)^. (19) 

To combine GRB data with the SNe la data and constrain 
cosmological models, we follow the simple method of avoiding 
any correlation between the SNe la data and the GRB data: the 
40 SNe points used in the interpolating procedure are excluded 
from the Union2 SNe la sample used to derive the joint con- 
straints (Liang, Wu and Zhang 2010; Liang, Wu and Zhu 2010). 
The 42 GRBs and the reduced 517 SNe la, CMB, BAO are all 
effectively independent, therefore we can combine the results by 
simply multiplying the likelihood functions. The total x 2 with 
the SNe + GRBs + CMB + BAO dataset is 

X 2 - /ffsNe.GRBs} + #CMB + AbAO' (20) 

We perform a global fitting to determine the cosmologi- 
cal parameters using the Markov Chain Monte Carlo (MCMC) 
method. In adopting the MCMC approach, we generate using 
Monte Carlo methods a chain of sample points distributed in 
the parameter space according to the posterior probability, us- 
ing the Metropolis-Hastings algorithm with uniform prior prob- 
ability distribution. In the parameter space formed by the con- 
straint cosmological parameters, a random set of initial values 
of the model parameters is chosen to calculate the^ 2 or the like- 
lihood. Whether the set of parameters can be accepted as an ef- 
fective Markov chain or not is determined by the Metropolis- 
Hastings algorithm. The accepted set not only forms a Markov 
chain, but also provides a starting point for the next process. We 
then repeat this process until the established convergence accu- 
racy can be satisfied. The convergence is tested by checking the 
so-called worst e-values [the variance(mean)/mean(variance) of 
1/2 chains] R - 1 < 0.005. 

Our MCMC code is based on the publicly available 
CosmoMC package (Lewis & Bridle 2002), and we gener- 
ated eight chains after setting R — 1 = 0.001 to guarantee the 
accuracy of this work. We show the 1-D probability distribu- 
tion of each parameter in the MCMC method (Q b A As, a, 
Qa, Age/Gyr, Q m , Hq) and 2-D plots for parameters between 
each other for the GCG model with SNe + GRBs + CMB 
+ BAO in figure 1 (Age/Gyr is the cosmic age, in units of 
Gyr). The best-fit values of the GCG model parameters with 
the joint observational data are A s =0.7475^|*(lo-)^™(2cr), 
a=-0.0256+° |™(lo-)+J} 2 ™(2cr), and the effective matter den- 
sity Q m = 0.2629+[5[5||(lo-)+°°2%(2cr). For comparison, fitting 
results from the joint data of 557 SNe la, the CMB and BAO 
without GRBs and 42 GRBs, the CMB and BAO without SNe la 
are given in Figs 2 and 3. We present the best-fit values of each 
parameter with the 1-cr and 2-cr uncertainties, as well as^ 2 . , in 
Table 1. 



From Figs. 1-3 and Table 1, it is shown that the cosmological 
constant (a = 0) is allowed at the 1-cr confidence level, and 
the original Chaplygin gas model (a — 1) is ruled out at 95.4% 
confidence level, which are both consistent with that obtained 
in Wu and Yu (2007c), and Li, Wu and Yu (2009). We can find 
that GRBs can provide strong constraints when combined with 
CMB and BAO data without SNe la, which has been also noted 
by Liang, Wu and Zhu (2010), Liang and Zhu (2010), and Gao et 
al. (2010). In additoin, the constraining results in this work with 
the joint observational data including GRBs are more stringent 
than previous results for constraining GCG model parameters 
with GRBs and/or other combined observations (e.g. Wang et al. 
2009a, 2009b; Freitasa et al. 2010; Davis et al. 2007; Wu and Yu 
2007c; Li, Wu and Yu 2009; Li, Li and Zhang 2010; Xu and Lu 
2010). 
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Fig. 1. The 2-D regions and 1-D marginalized distribution 
with the 1-cr and 2-cr contours of parameters Q b /i 2 , As, a, 
Qa, Age/Gyr, Q m , and Hq in GCG model, for the data sets 
SNe+GRBs+CMB+BAO. 
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Fig. 2. The 2-D regions and 1-D marginalized distribution 
with the 1-cr and 2-cr contours of parameters Q.\,h 2 , As, a, 
Qa, Age/Gyr, Q m , and Hq in GCG model, for the data sets 
SNe+CMB+BAO. 
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The GCG Model 
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Table 1. The best-fit values of parameters Q,\>h 2 , As, a, Qa, Age/Gyr, Q,„, and Hq for the GCG model with the 1-cr and 2-cr 
uncertainties, as well as^ in , for the data sets SNe+CMB+BAO, SNe+GRBs+CMB+BAO, and GRBs+CMB+BAO, respectively. 
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Fig. 3. The 2-D regions and 1-D marginalized distribution 
with the 1-cr and 2-cr contours of parameters £l\,h 2 , As, a, 
Qa, Age/Gyr, Q m , and Hq in GCG model, for the data sets 
GRBs+CMB+BAO. 



4. CONCLUSIONS 

By using the Markov Chain Monte Carlo method, we have 
constrained on the generalized Chaplygin gas (GCG) model 
with the cosmology-independent GRBs, as well as the Union2 
SNe la set, the CMB observation from WMAP7 result, 
and the BAO observation from SDSS DR7 galaxy sample. 
With the joint observational data, the best-fit values of the 
GCG model parameters are A s =0.7475!^|*(lo-)!^^(2cr), 
a=-0.0256+° |™°(lo-)+Jj-2730 (2o .^ and the effective matter den- 
sity Q m = 0.2629+°-°!|(l(x)+°-°^(2cr), which are more strin- 
gent than previous results for constraining the GCG model pa- 
rameters obtained using data of GRBs and/or other combinations 
of observations. 
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